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Xyy I Abstract. A new scheme of radiation transfer for understand- 

^—^ • ing infrared spectra of H II regions, has been developed. This 

CN i scheme considers non-equilibrium processes (e.g. transient heat- 

^^ ' ing of the very small grains, VSG; and the polycyclic aromatic 

^ . hydrocarbon, PAH) also, in addition to the equilibrium thermal 

C^^ I emission from normal dust grains (BG). The spherically sym- 

''^ ' metric interstellar dust cloud is segmented into a large number 

Q>^ . of "onion skin" shells in order to implement the non-equilibrium 

^P I processes. The scheme attempts to fit the observed SED origi- 

nating from the dust component, by exploring the following pa- 
rameters : (i) geometrical details of the dust cloud, (ii) PAH size 
and abundance, (iii) composition of normal grains (BG), (iv) 
radial distribution of all dust (BG, VSG & PAH). 

The scheme has been applied to a set of five compact H II re- 
^ ; gions (IRAS 18116-1646, 18162-2048, 19442+2427, 22308+5812, 

Cd ■ & 18434-0242) whose spectra are available with adequate spec- 

* * I 

^ , tral resolution. The best fit models and inferences about the 

l>( I parameters for these sources are presented. 

C^ ' Key words: H II regions - radiative transfer - PAH - VSG 



1. Introduction 

Till recently, the mid to far infrared spectral energy distribution (SED) of 
Galactic star forming regions in general was available only in the four IRAS 
bands (12, 25, 60 &: 100 fim). In some relatively rare cases, spectroscopy 
in the 10 fim. band through the atmospheric window, was also available. 
However, the situation has changed drastically recently, due to the advent 
of the Infrared Space Observatory (ISO). The ISOPHOT photometer along 
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with ISO-SWS and ISO-LWS spectrometers together has revolutionized the 
availabihty of information about SED of the astrophysical sources in general. 

In the literature, several radiation transfer schemes have been used for 
interstellar dust clouds with embedded YSOs in spherical (e.g. Scoville &: 
Kwan 1976, Leung 1976, Churchwell, Wolfire, & Wood 1990) as well as cylin- 
drical (Ghosh & Tandon 1985, Dent 1988, Karnik & Ghosh 1999) geometries. 
All of these considered the dust grains to be in thermal equilibrium. The 
role of non-equilibrium processes (resulting in transient heating / excitation 
of grains, particularly in the vicinity of a source of UV radiation) has become 
evident from significant near & mid infrared continuum emission detected 
in Galactic star forming regions (Sellgren 1984, Puget, Leger, & Boulanger 
1985, Boulanger, Baud & van Albada 1985) as well as spectral features 
(Leger &: Puget 1984, Allamandola, Tielens, & Barker, 1985, Puget & Leger 
1989). The importance of these processes has also been demonstrated in ex- 
tragalactic nuclei / star forming regions (Moorwood et al., 1996, Metcalfe et 
ai, 1996, Ghosh, Drapatz, & Peppel, 1986). A comparison of the observed 
mid-IR spectral features with laboratory data, has led to the identification 
of a new constituent of the interstellar medium - polycyclic aromatic hydro- 
carbons (PAH). The enhanced continuum emission in the near & mid IR has 
been mainly attributed to the very small grains (VSG) of radii 10-100 A. 
Hence, it is obviously important to include the non-equilibrium processes in 
attempting to model the observed SED of star forming regions in general. 
Basically, grains of very small size or a large organic molecule, with effective 
heat capacity comparable to the energy of a single UV photon get excited 
(for a short time) to an energy state well above its thermal equilibrium 
state corresponding to the local radiation field. The photons emitted during 
the de-excitation process contribute to the near / mid IR part of the SED, 
which shows continuum excess and emission features which are unexplained 
by radiative transfer models considering the emission from large grains in 
thermal equilibrium alone. Recently, Siebenmorgen & Krugel (1992) have 
attempted to quantify the properties of the dust components relevant for 
non-equilibrium processes (VSG & PAH), from the infrared data of sources 
in different astronomical environments in our Galaxy. The role of VSG on 
the infrared emission from externally heated dust clouds has been studied 
by Lis & Leung (1991). Krugel &: Siebenmorgen (1994) have presented a 
method to model the transfer of radiation in dusty galactic nuclei, which 
includes the presence of VSG and PAH. 

Here we present a scheme of radiative transfer developed by us which 
is applicable in spherical geometry. This includes, in addition to the dust 
grains in thermal equilibrium (of normal size, hereafter big grains or BG), the 
transient heating of very small grains (VSG) as well as the PAH molecules. 
An attempt has been made to model five compact H II regions : IRAS 
18116-1646, 18162-2048, 19442+2427, 22308+5812 & 18434-0242 using the 
above scheme. 
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In section 2, the radiative transfer modelling scheme is briefly described. 
The results of modelling the five compact H II regions are presented in 
section 3. The last section (4) consists of discussion. 



2. The Modelling Scheme 

2.1 Dust components and their properties 

The normal grains (BG) consist of two components : astronomical silicate 
and graphite. Their size distribution is taken as per Mathis, Rumpl &: 
Nordsieck (1977) to be a power law, n{a) oc a'^, with -3.5 as the exponent. 
The lower and upper limits of the grain radii are taken to be 0.01 ^m and 
0.25 ^m respectively as recommended by Mathis, Mezger &: Panagia (1983), 
for both astronomical silicate as well as graphite grains. The scattering and 
absorption coefficients, and anisotropic scattering factors have been taken 
from Draine & Lee (1984) and Laor & Draine (1993). 

The VSG component is taken to be graphite grains of a single size : 
either 10 A or 50 A in radius. Their optical properties have also been taken 
from Draine & Lee (1984). Abundance of VSG is connected to that of the 
normal grains through a scaling factor YvsGi which gives the fraction of dust 
mass in VSG form to the normal BG form. The value of YvsG was taken 
from Desert, Boulanger, & Puget (1990), which is needed to account for the 
2200 A bump in the average interstellar medium in the Galaxy, and it has 
been held fixed for all models considered here. 

The PAH component is assumed to be either a single molecule with about 
15 - 30 atoms, or a large complex consisting of 10 - 20 of these molecules as 
used by Siebenmorgen (1993). Their optical properties, feature centres, fea- 
ture shape and widths have been taken from Leger & d'Hendecourt (1987). 
The abundance of PAH component is also connected through a scaling fac- 
tor YpA_H to the normal grains (BG). There are two additional parameters 
explored in the modelling of the observed PAH spectral features : (i) the 
radius of the PAH molecule / complex, a pah'-, and (ii) the de-hydrogenation 
factor, fde-H- The value of fde-H lies between and 1 {fde-H = refers 
to completely hydrogenated PAH). Whereas apAH has implications of heat 
capacity and hence the efficiency of transient heating for a given radiation 
field, the fde-H affects the ratios of PAH features resulting from the C-H 
versus C=C stretch modes. 



2.2 Geometry 

The star forming region is considered as a spherical dust cloud immersed in 
an isotropic interstellar radiation field, with an embedded source of energy 
(e.g. a ZAMS star) at its centre. A central cavity in this cloud represents 
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Figure 1. Schematic diagram of the shell structure of the cloud. 



sublimation / destruction of grains in the intense radiation field of the central 
source. A schematic of the dust cloud is presented in Figure 1. 

This spherically symmetric dust cloud, is divided into a large num- 
ber of concentric contiguous spherical shells (say Shi, Sh2, ■■■, Sh^) hke 
"onion skins". Each shell, Shi, is identified by its inner and outer radii 
(^j^rnm ^ RV^"-^ ; gee Figure 1). These sheUs can be of different selectable 
thicknesses, depending on the optical depth at the shortest relevant wave- 
length. In order to incorporate the presence of both - normal grains (BG, 
responsible for emission at thermal equilibrium), as well as the grains respon- 
sible for non-equilibrium emission (VSG and PAH), each shell is subdivided 
into a pair of sub-shells, Shf^ & ShJ^ corresponding to these two com- 
ponents respectively. Whereas the former consists of only BG, the latter 
consists of only the VSG and PAH. 

The full detailed radiative transfer calculations assuming the normal 
grains to be in thermal equilibrium, are performed in each of the sub-shells 
Shf , for i = 1,2,..., N. The subshells ShJ^, go through a statistical me- 
chanical treatment describing the non-equilibrium emission processes for the 
VSGs and the PAHs. For simplicity of computations, the sub-shells Shf 
are considered to be very thin compared to the total thickness of the shell 
Shi ) and this sub-shell is assumed to be placed at the inner edge of the shell 
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Shi (see Figure 1). The final results are expected to be insensitive to the 
above simplification since individual shells are optically thin. 

Radiative transport at each of the two sub-shells is carried out as a 
two point boundary value problem, the two boundary conditions being the 
incident radiation fields at the two surfaces. The calculations begin with 
the given spectrum emitted by the embedded energy source (in general, an 
Initial Mass Function weighted synthetic stellar spectrum ensemble) incident 
at the inner boundary of the first shell Shi. The outer surface of the last 
(outermost) shell, Sh]^, has the interstellar radiation field (ISRF) incident 
on it from the outside. Starting from the "core" side of the first shell, the 
radiation is transported through the sub-shell ShY^ first and the emergent 
processed spectrum is considered to be incident on the other sub-shell Shf^ . 
The emergent spectrum from the latter is the processed output of the entire 
shell Shi and is used as input boundary condition for the next shell Sh2- In 
this manner, the radiation field is transported outward from shell Shi to Sh2 
... till the last shell, viz., Shjsf is reached. This entire processing from shell 
1, 2, ... to N, constitutes one iteration. Several such iterations (typically 
5-10) are carried out until a set of predetermined convergence criteria are 
satisfied. The emerging spectrum from the last shell, Sh]\f, is the desired 
output of the full model. The number of shells used for a specific source is 
determined by the criterion that the shell is optically thin in the shortest 
relevant wavelength. 

2.3 Processing of transient heating of the VSG and the PAH 

As described above, the dust components in a shell for which the non- 
equilibrium emission processes are important, are segregated into a sepa- 
rate sub-shell (Sh^^) consisting only of the VSG and the PAH. The in- 
teraction of the total incident radiation field from both the surfaces of this 
sub-shell, {Iin)i with the VSG component, is considered using a code de- 
veloped by us based on the statistical mechanical treatment prescribed by 
Desert, Boulanger & Shore (1986). The incident radiation, partly extin- 
guished by the VSGs {lyp = I in x e^'^v'sc), is considered incident on the 
PAH component and a similar computation is repeated. The final emerging 
spectrum consists of three components: (i) the originally incident radiation 
extinguished by both VSG as weU as PAH, {I^ut = lin x c'^^^^^sg+^pah)), 
(ii) the emission from the VSG component, and (iii) the emission from the 
PAH component. Whereas the first component is direction sensitive (the 
two surfaces get different contributions depending on the original spectrum 
incident at the other surfaces), the latter two contribute equally to the two 
surfaces. 

The VSG and PAH components of grains have fiuctuating temperature, 
mainly because their enthalpy (internal energy) is comparable to the energy 
of UV or visual photons. This means, the multiphoton absorption processes 
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can become important (depending on the exact radiation field and the de- 
tails of thermal & optical properties of these grains) as they can lead to a 
modified temperature distribution. An iterative method has been used here 
to consider these multiphoton processes for VSGs and PAHs separately. The 
method assumes a single grain in an isotropic radiation field, and follows the 
evolution of the grain temperature by solving the relevant stochastic differ- 
ential equation. 

A scheme of between 100 to 400 levels of internal energy (covering 0.5 
eV to 200 eV) for considering discrete heating / cooling processes; and 400 
energy levels (for energies 1.25 xlO~^ eV to 0.5 eV) for considering the 
continuum processes, has been incorporated. A total of 97 frequency grid 
points covering 0.0944 /Ltm to 5000 ;um have been used. Several grid points 
are densely packed around the five PAH features at 3.3, 6.2, 7.7, 8.6 and 
11.3 //m. 

2.4 Radiation transport through normal grains (BG) 

Each of the sub-shells consisting of the normal grains, {Sh^ , k = 1,2, ...N), 
separately undergoes full radiative transport calculation using the code CS- 
DUST3 developed by Egan, Leung & Spagna (1988) (see also Leung 1975). 
In CSDUST3, the moment equation of radiation transport and the equation 
of energy balance are solved simultaneously as a two-point boundary value 
problem. The effects of multiple scattering, absorption and re-emission of 
photons on the temperature of dust grains and the internal radiation field 
have been considered self-consistently. In addition, multi grain components, 
radiation field anisotropy and linear anisotropic scattering are also incorpo- 
rated. 

Same frequency grid of 97 points, as used for VSG and PAH, has been 
used here. In order to avoid non-convergence problems due to sharp changes 
in optical depth at any of the frequency grids, logarithmically increasing 
radial grid spacings have been used at the inner shell boundary. Similarly 
a smoothly decreasing grid spacings have been used near the outer shell 
boundary. 

2.5 The Modelling Scheme 

The scheme aims to construct a model constrained by the observed SED 
covering the entire infrared and the sub-mm / mm region. Based on com- 
parisons of the model predicted SEDs with the observed SED, various model 
parameters are fine tuned till the best fit model is identified. The following 
model parameters are explored : (i) the total radial optical depth (repre- 
sented at a fiducial wavelength of 100 /im); (ii) exponent of the dust density 
distribution power law; (iii) the ratio of graphite component to the astro- 
nomical silicate component for BGs; (iv) size of VSGs, aysG'- either 50 A or 



Radiation transfer in H II regions 



67 



Table 1. Input parameters of the compaet H II regions 



IRAS Source 


L 


T, 


-L'sun 


Odia 


Rout 




(Lg) 


(K) 


(kpc) 


n 


(pc) 


18116-1646 


1.6 xlO^ 


40,000 


4.4 


105 


1.12 


18162-2048 


2.8 xlO^ 


30,900 


1.9 


63 


0.29 


19442+2427 


5.4 xlO-^ 


35,500 


2.3 


93 


0.52 


22308+5812 


8.8 xlO^ 


37,500 


5.7 


95 


1.31 


18434-0242 


1.0 xlO^ 


48,000 


7.4 


38 


0.69 



10 A ; (v) size of PAH molecule / cluster, qpah'- 4.6 A or 8 A or 13.6 A ; (vi) 
relative abundance of PAH compared to BGs, Ypah (a constant or varying 
with the radial distance) ; and (v) the de-hydrogenation factor, fde-H- 

The inner radius (Rin), for each model of the spherical cloud, has been 
determined using the constraint that the temperature of the BG is equal 
to 1500 K, the sublimation temperature of the normal big grains (graphite 
and astronomical silicate). The radial dust density distribution law has been 
assumed to be a power law Udir) oc r^", and the values of a that have been 
explored are 0, 1 and 2. 

3. Application of the modelling scheme 

In order to demonstrate the usefulness of the scheme described above, an 
attempt has been made to apply the same to a few H H regions. The 
question is : in spite of rather simplistic treatment, can we get any insight 
into physical details of these sources? 



3.1 The sample of compact H II regions 

With the advent of Infrared Space Observatory (ISO), it has become now 
possible to have precise spectroscopic information in the entire infrared band 
encompassing near to far infrared region. The spectroscopic results for a 
sample of six Galactic compact H II regions, covering four of the five major 
PAH features have been published by Roelfsema et al. (1996). We have 
chosen five out of their six sources for our detailed study. The sixth source 
IRAS 21190+5140, identified with Ml-78 and variously considered as H II 
region and planetary nebula (Puche et al. 1988, Acker et al. 1992), has not 
been considered here. Although the published spectral results from ISO is 
rather limited (6 - 12 //m), if the IRAS Point Source Catalog (IRAS PSC) 
measurements (at 12, 25, 60 & 100 jum), IRAS Low Resolution Spectra 
(IRAS LRS; between 8-22 /xm) and ground based spectroscopy around 
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Table 2. Dust parameters valid for the entire sample of compact H II regions 



Dust component 


Parameter (unit) 


Value 


BG: 


arnin (^m) 


0.01 




ttrnax (^m) 


0.25 




7 


-3.5 


VSG: 


avsG (A ) 


50.0 




YvSG 


4.70 xlO-^ 


PAH: 


apAH (A ) 


8.0 




YPAH 


4.30 xlO-3 




fde-H 


0.0 



the 3.3 /um PAH feature, are included (whenever available), then sufficient 
observational constraints can be placed on the radiative transfer models. All 
these measurements have been compiled to construct the SEDs for the five 
compact H H regions, which are displayed in Figure 2. 

The total luminosity, Lfot, has been taken from Roelfsema et al. (1996) 
for all the sources except for IRAS 18434-0242. The luminosity for IRAS 
18434-0242 listed by them is too low compared to that estimated from the 
IRAS data. The latter has been used here for modelling. Assuming the 
embedded source to be a single ZAMS star of luminosity Ltot, a Planck- 
ian spectral shape with corresponding temperature taken from Thompson 
(1984), has been assumed. Since distance estimates are available for all these 
sources, angular sizes are required to fix the outer radii of the clouds. The 
mid infrared angular sizes are estimated by comparing the flux densities at 
12 fu.ni as measured by IRAS-LRS and the ISO-SWS and the solid angle 
covered by the ISO-SWS. It is assumed that the source size is smaller than 
the IRAS-LRS beam and has constant brightness per unit solid angle. The 
entire size of the cloud has been estimated from the 12 ^m size by using 
empirical relation between angular size and the wavelength for compact H 
II regions (Mookerjea & Ghosh 1999). 



3.2 Results of modelling 

It has been possible to get reasonable fits to the observed SEDs of all the 
five compact H II regions by varying parameters of our modelling scheme. 
The predicted spectra and the observations are compared in Figure 2. 

The following comments are valid for all the five compact H II regions 
studied here. It was found that the models with uniform density distribu- 
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Figure 2. Spectral energy distribution of the compact H II region IRAS 18116- 
1646. The ordinate is the log of the flux density multiplied by the surface area 
of the respective cloud. Solid lines show the ISO-SWS spectra from Roelfsema 
et al. (1996); dotted lines show our best fit model spectra; diamonds show other 
observations. Other observations include — 3 /im observations from de Muizon, 
d'Hendecourt & Geballe (1990), IRAS LRS spectra from Olnon & Raimond (1986) 
or from Volk & Cohen (1989), IRAS PSC flux densities, sub- mm observations from 
McCutcheon et al. (1995), Jenness, Scott, & Padman (1995) or Barsony (1989), and 
1.3 mm observation of Chini, Krugel, & Kreysa (1986). In order to show the PAH 
features clearly, mid IR region of the SEDs are shown separately (on the right). 



tion i.e. n{r) oc r" (as opposed to n{r) oc r^^ or r^^) gave much better 
fits to the SEDs. The VSGs with aysG = 50 A and the PAHs with inter- 
mediate size (i.e. apAH = 8 A ) give better fits to the respective spectra. 
The de-hydrogenation factor, fde-H-, needs to be zero (corresponding to 
Nh = \/6 X Nc) in order to fit the relative strengths of the PAH features 
for all the five sources. This value of fde-H-, is typical for the types of PAH 
which have been strongly proposed in the literature, viz., Coronene & Ova- 
lene (Leger & Puget, 1984). In addition, this fde-H is consistent with the 
value of apAH (8 A ) inferred from our modelling, since such PAHs are ex- 
pected to be completely hydrogenated in the emission zones (Allamandola, 
Tielens & Barkar, 1989). Table 2 lists those best fit parameters related to the 
various dust components, which are valid for the entire sample of compact 
H H regions. 

It has been found that, whereas the BG and VSG components should ex- 
ist throughout the cloud, it is absolutely necessary that the PAH component 
must be confined to a thin inner region {Rg^^^ << Rout), in order to repro- 
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Figure 2. Continued ... (for IRAS 18162-2048) 
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Figure 2. Continued ... (for IRAS 19442+2427) 



duce the PAH features. The size of this region is quantified by a parameter 
'nPAH-, which is defined as : 'qpAH = {{Rot^^ - Rin)/{Rout - Rin))- This 
parameter had to be varied for each source, till a good fit to the spectrum 
was obtained. In addition, the abundance of PAH relative to BGs, Ypah^ 
needed to be increased by a factor 10 relative to the normal value obtained 
by Desert, Boulanger & Puget (1990). However this does not lead to any 



Radiation transfer in H II regions 



71 









Fig 


2 (4_a) 




26 








^,,--f"' 


'22308+5812 MODEL' 

■22308+5812 SWS' 

'22308+5812_OTHERS' + " 


25 






j/' 




■■■■., 


^ 24 

Q 
m 23 






f 




\, 


f 


- i^ 










22 








■- 


21 










- 



'22308+581 2_MODEL' 

22308+581 2_SWS' 

■223D8+5812_OTHERS' 




>V^^"^'^ 



Figure 2. Continued ... (for IRAS 22308+5812) 



.i^ 



'18434-0242_MODEL' 
'18434-0242_SWS' ■ 
'18434-0242_OTHERS' 



log-lamda(micron] 



Figure 2. Continued ... (for IRAS 18434-0242) 
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conflict with the available carbon, since t/pah << 1. The values of the best 
fit parameters specific to each source, are presented in Table 3. 
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Table 3. Best fit parameters of tfie compact H II regions as determined by mod- 
eliing 



IRAS Source Name 


Rin 


Rout 


nn 


MtoI 


^tot 
^100 


Graphite : Silicate 


VPAH 




(pc) 


(pc) 


{cm-^) 


(Mo) 




(% : %) 




18116-1646 


1.4 xlO-3 


1.12 


1.32 X 10^ 


1.9 X 10^ 


0.056 


75 : 25 


2.1 xlO-2 


18162-2048 


5.3 xlO-4 


0.29 


1.32 X 10^ 


3.3 X 102 


0.14 


88 : 12 


2.7 xlO-2 


19442+2427 


7.5 xlO-4 


0.52 


5.30 X 10^ 


7.7 X 102 


0.10 


95 : 5 


1.5 xlO-2 


22308+5812 


1.5 xlO"-'^ 


1.31 


1.38 X 10*^ 


3.2 X 10^ 


0.068 


77 : 23 


1.3 xl0"2 


18434-0242 


2.2 xl0"3 


0.69 


5.30 X 10^ 


1.8 X 10^ 


0.14 


95 : 5 


5.5 xlO-2 



Discussion 



The following inferences can be drawn about the sources modelled here, 
provided the basic assumptions (e.g. spherical symmetry; sources of energy 
located only at the centre of the cloud; etc) are not at great variance from 
the reality. 

The most favoured radial dust density distribution law, for all five sources, 
turns out to be of uniform density. This can perhaps be understood in terms 
of the far infrared constraints (IRAS-PSC 60 & 100 /um data). If the dust 
density is falling with radial distance, then in order to fit the FIR part of 
the SED, so high a dust density is required at the vicinity of the embedded 
ZAMS star, that the mid infrared emission becomes invisible. This prob- 
lem can perhaps be avoided in a non-spherically symmetric geometry. We 
have explored the effects of relaxing the assumption that Rin = Ro, where 
Rq is the radial distance at which (BG) grain temperature becomes equal to 
the sublimation temperature (1500 K). The most important effect of making 
Rin > Rq, is to drastically modify the near and mid infrared continuum level 
of the predicted spectrum. In addition, the role of non-equilibrium processes 
vis-a-vis thermal equilibrium emission of the PAH features, as well as the 
continuum due to VSG, changes significantly. 

A quick perusal of Figure 2 and Tables 2 and 3 brings out the following 
facts: 

1) All the compact H II regions considered here, are deeply embedded 
stars; total optical depth at 100 /im in the range of 0.056 - 0.14. This is 
necessary to explain the far IR spectra observed by IRAS. 

2) PAH is confined only to a thin central shell; the thickness of this 
shell being just a few percent (1.3 - 5.5 %) of the total thickness of the 
dust cloud. As these sources are optically thick at mid IR, if the PAH is 
distributed throughout the cloud, its emission which occurs in the inner hot 
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region where high energy photons responsible for non-equihbrium processes 
are present, will be absorbed by the outer cooler shells, and PAH features 
will not be detectable. 

3) The BGs are dominated by graphites, with silicates contributing less 
than 25 %. The latter has been tied down rather precisely by the 10 fiui 
silicate feature. 

4) ISO-SWS fluxes are generally much smaller than IRAS fluxes at similar 
wavelengths, indicating that SWS is not sampling full emission at mid IR and 
the source sizes at these wavelengths are much larger than the SWS beam 
size (14" X 20"). With this in mind, we have not tried to fit the absolute 
fluxes of the SWS but only used its shape as indicative of the importance of 
PAH molecules. 

Following comments can be made about the individual sources: 

• IRAS 18116-1646 : It has relatively lower optical depth. The fit to 
the IR data is quite reasonable, except at 100 ^m where IRAS flux is 
higher; no sub-mm observation exists for this source. 

• IRAS 18162-2048 : This source (GGD27) was originally thought 
to be a HH object. However now it has been established as a star 
forming region with reflection nebulosity as well as outflow (see for 
example Stecklum et al. 1997). The region has several near IR and 
mid IR sources; the source of energy being close to IRS2. The size of 
this source at sub-mm wavelengths is ~ 1' (McCutcheon et al. 1995), 
consistent with the size for the best fit model. The mass of the envelope 
estimated by our model is not far from the estimate of Yamashita et 
al. (1987), viz., 2OOM0. They have proposed a disk geometry for this 
source. This source has very high optical depth. The fit to the IR data 
is quite reasonable but at the sub-mm wavelengths the calculated flux 
densities are lower than the observed ones. 

• IRAS 19442-1-2427 : This source lies in the H II region S87. The 
size of this source at sub-mm wavelengths is ~ 1' (Jenness, Scott and 
Padman, 1995), consistent with the size for the best fit model. This 
source has medium optical depth. The fit to all the observations from 
3 /im to 850 ;um is quite reasonable. 

• IRAS 22308-1-5812 : It has relatively lower optical depth. The fit to 
the IR data is quite reasonable; no sub-mm observation exists for this 
source. 

• IRAS 18434—0242 : This source is the most luminous source with 
high optical depth. There are no IR observations for this source other 
than those from IRAS and ISO. IRAS PSC 100 /im as well as 1.3 mm 
observations are higher than calculated. 



74 S.K. Ghosh & R.P. Verma 

From the above, we conclude that our new scheme of radiative transfer 
which includes non-equlibrium processes (transient heating of the grains / 
PAH / VSG) in addition to the emission in thermal equilibrium, can give 
important physical insight into Galactic star forming regions. If the simplify- 
ing geometrical assumptions of our scheme are valid, then several important 
inferences can be made about the five compact H II regions considered for 
modelling here. 
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